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ON SINGULAR POINTS IN THE APPROXIMATE DEVELOPMENT 
OF THE PERTURBATIVE FUNCTION* 

BY 

FRANK WALKER REED 

§ 1. Introduction. 

A method of computing approximately the coefficients of those terms of high 
order in the perturbative function which derive their importance from the pres- 
ence in them of a small divisor has been set forth by Poincar6. As applica- 
tions of the method to the problem of three bodies, Poincar:e; f has considered 
the case of coplanar orbits, one elliptic and one circular, and the case of 
coplanar elliptic orbits, both eccentricities being small. The latter case has 
been developed in detail by Coculesco. J Hamy § has discussed the cases of 
an elliptic enclosing a circular orbit and a circular enclosing an elliptic orbit, the 
mutual inclination being small, and the case of circular orbits with any inclina- 
tion. The latter of these was discussed anew by F:fiRAXJD || together with the 
case of an elliptic enclosing a coplanar circular orbit, and the case of a circular 
enclosing a coplanar elliptic orbit. 

Owing to the high degree of the algebraic equations which define the singu- 
lar points upon which the approximate values of the coefficients depend, it is 
found necessary, in the general case, to resort to numerical examples. This 
paper is devoted to a discussion of these algebraic equations in certain cases 
more general than those so far considered, and to an extension of the criterion 
of admissibility as employed by PoiNCAR^. The singular points in the most 
general case and in an important particular case are found and classified. The 
results in the cases studied by P01NCAR6 and Coculesco are given first because 
of their usefulness in treating the more general problem. Before proceeding to 

* Presented to the Society (Princeton) September 13, 1909. 

t PoiNCAEB. Les Methodes Nouvelles de la Mecanique CSleste, vol. 1, Chap. VI. 

t Coculesco. Sur les expressions approchees des termes d'ordre eleve dans le dSveloppemenl de la 
fonciion perturbairice, Journal de Mathematiqnes, ser. 4, vol. 1 (1895). 

§ Hamy. Sur les ineqalites d'ordre 6leve, Bulletin Astronomique, vol. 10 (1893). Sur le 
developpemeni approche de la fonciion perturbatrice, Journal de Mathematiques, ser. 4, vol. 
10 (1894) ; ser. 5, vol. 2 (1896). 

II Feeaud. Sur la valeur approche des coefficients d'ordre Ueve dans le developpemeni en series, 
Annales de I'Observatoire de Bordeaux, vol. 7. 
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486 F. w. reed: on singular points in the approximate [October 

this discussion, it is necessary to set forth briefly the principles of Poincar6's 
method, but in doing so the general equations are derived at once. 

It is required to develop a certain function F\ in terms of the mean anom- 
alies I and I' in the form 



00 00 



^? = E E ^M,M,e"-''''^'^''\ i=V-l, and m,, m, integers. 

mi—Q m2=:0 

The coefficients A^^^^ are given by the integral 

where 
and 

The new variables, t and z, have been introduced by the transformation 

e*^ = t% e*'' = t-''z^, 

and the new constants by the relations 

m^ = ati + h, m^ = en + d, 

where a, b, c, d are finite integers, small relative to w, and a and c are rela- 
tively prime. Now if $(«) is expanded as a Laurent series 

and the result substituted in the expression defining A^^^, we get 

A — n 

A theorem by Darboux * can be employed to obtain approximate values of a^ , 
and consequently of A^^^ , for large values of n without computing the interven- 
ing coefficients. 

The initial contour of integration, | ^ | = 1 , may be deformed in any manner 
without affecting the value of the integral <E> (s ) provided no singular point of 
F{z, t), considered as a function of t with s as a parameter, is encountered. 
By varying z continuously the singular points of F(z, t) are displaced continu- 
ously. Singular points of *!*(») are those values of z for which two singular 

*Daeboux. Sur I' approximation des fonctions de ir^s-grands nombrea, Journal de Math4- 
matiqaes, ser. 3, vol. 4 (1878). 
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points of F{z, i), considered as a function of t, come into coincidence. If the 
points were initially, that is for | a | = 1 , one inside and the other outside the 
contour |^| = 1, then the singular point a of ^{z) is admissible ; for, ^(») 
ceasing to be holomorphic, the theorem of Darboux is applicable. 
The principal part of the perturbative function * is 

where 



2 



yS = i''i-^e'("'-">, ^, = L''L 



-m) 



I/, aia (j), I, u, a> represent respectively the major axis, the eccentricity, the 
mean anomaly, the eccentric anomaly, and the longitude of perihelion of a planet 
B referred to the line of nodes and the center of gravity, D, oi A and i? in a 
system of three bodies A, B, C. L'^, sin <f>', etc. are the elements of C 
referred to D and the lines of nodes, t/is the inclination. ^, rj are defined by 

^ = cos u — sin (f> + i cos (j) sin u , 

t] ^ cos u — sin 4>' + i cos (j)' sin u, 

and f^, rjg&ve their conjugates. With this notation the coordinates oi B and C 
respectively are given by the real and imaginary parts of 

respectively. 

The singular points of i''(s, t) with which we are concerned are the same as 
those of i^J. These are (a) the points for which u and u' as functions of I and 
r cease to be uniform, given by 

dl . 5Z' ^ . , , ^ 

-^^r— ^ 1 — sin rf) cos u ^ 0, -^— =1 — sin 4> cos u ^ , 

du ^ ou 

and (6) the singularities of F\, given by 

P' = 0. 
Making the substitutions 

x = e'», y=e"", i = |(cosJ'-l), 
tan ^- = T, 1 + T^^tr, 

tan ^ = T, l+T =o", 



* This form of the perturbative function is due to Feeaud, loo. oit. 
Trans. Am. Math. Soc. 38 
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we get 

{x-rl (1-xr)^ 

^~ o-tc ' ^»~ o-a; ' 

{y-r'f {i-yr'f 

'^- T'y ' ""- T'y • 
Then the equations defining the singularities become 

(1) (cc-t)(1-xt) = 0, 

(2) (2/-t')(1-2/t') = 0, 
P = {aa'xyyr = \_(T'y{x - rf - ^<Tx{y - r'f ] 

(3) X \(T'y{l—XTf — ^^(Tx{l—yr')-'\ + ja-a'xy[y{x — rf(y — r'f 

+ 7o(l- xt)^(1- 2/t7- ^„(x - T)X1- 2/t7- ^(1 - XT)Xy - r'f] = 0. 

Singular points of the Jlrst species of ^{z) are given by the intersections of 
the curves (1), (2), (3). Singular points of the second species are those satisfy- 
ing the condition for double roots of (1), (2), or (3), any one of which may be 
written in the general form in the original variables 

f{z,t) = 0. 

The condition for a multiple root ofy(s, t) is 

•^ ' dt {x — r) {1 — xr) dx (y — t')(1 — yr' ) dy 

Equations (1) and (2) will not give rise to points of the second species. 
Consider equation (3), which may be written 

■^(^K' y) = BC+jaa'xyA, 
where 

A=^j{x- rf{y -r'f + 7„(1 - xry(l - yr'f 

-^^(x-r)\l-yr'y-^{l-xr){y-r'y, 
B = (T'y{x — ry — ^(7x{y — t')^, 
G = <7'y{l-ayry-^^ax{l-yr'y. 
The equation d P jdt = becomes, after reducing by means of P = , 
R = B' C{x - r){y - t) - B C {1 - xr){\ - yr) 
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-jc<7x{y-T'){l-yr')^A-x-^^ 

+ja(T'y{x-T){l-XT)\A-y-^^ =0, 

where 

B'=c{x + t)(1 - yr') + a^{y + t')(1 - xt), 

C = c{l + XT){y - r') + a^^{l + yT'){x - r). 

§ 2. CoPLANAE Orbits. 

The equations P =0 and H = 0, which determine all the singular points of 
the second species in the most general case, admit of important simplifications 
in certain particular cases. If we suppose the inclination zero, then j =0, and 
we have 

(1) (a!-T)(l-a5T) = 0, (3) B=0, (6) B' = 0, 

(2) (y-T')(l-2/T')=0, (4) (7=0, (6) C"=0, 

where (3), (4) are the factors of the reduced equation P = , and (5), (6), 
those of ^ = . In this case the singularities are as follows : 

B=0, x=T, and y=T' intersect in U, C=0 , and y=T' intersect in V and W\ 
(7=0, a;=-, and 2/=-, intersect in V", i?=0, and 3/=- intersect in F' and TF, 

(7= 0, and 35= T intersect in jE' and JT', »5=t, and 3/=— intersect in JT, 

5=0, and 05=- intersect in ^'and Z', a;=-andv=T' intersect in Jf', 

B = and (7=0 intersect in B', B, B, and B', 
B=0 and B=0 intersect in C", I), T, and F, 
(7= and (7'= intersect in T', B', C, and F'. 
These singular points have been discussed by Poincar^; and Coculesco. 

§ 3. The General Case. 

The main difficulty which arises in the general case is that of solving the 
algebraic equations. Even when t, t', j are small, as is usually the case in 
nature, the approximations thus made possible do not so reduce the equations 
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that they admit of literal solutions. But the problem is not without interest 
from the theoretical standpoint, although the method it is necessary to employ is 
dependent upon numerical assumptions. In fact the results to be obtained in 
this manner have the same generality as those which would be furnished by a 
literal solution. 

Since for j' = the equations P = and i? = reduce to 

5 = 0, (7=0, 5 = 0, C" = 0, 

it is easily seen that all the singular points for a small inclination will lie near 
those given in § 2. The reducible curve 5 = will lie near the curves 5 = 
and C = , but 5 = 0, being irreducible, will not lie near 5' = or C = . 
Referring to equations (3), (4) of § 1, and recalling the fact that /3, /S^ and 
also 7, 7j are conjugate complex quantities, we can easily verify that 

A{x,y,i) = x'y'A(^-, -, -. j, 
A — x-i^r- = A, (x, y, i) = — x^y^A, ( -, -, v | , 

^-y^y=^^^''^y^'^=-''y^\x^y^j)- 

Substituting in P and 5 above Ijx, l/y, 1/i for x,y,i, multiplying through- 
out by x*y*, and changing signs in the latter function, we arrive at the identical 
expressions for 5 and 5. Thus 5=0 and 5 = are self -reciprocal equa^ 
tions in x, y, provided that /3 , yS^ and also 7 , 7„ interchange. This was to 
be expected, since the equations were originally in cos u , cos u, sin u , sin u 
which involve x, y reciprocally. Consequently the discussion of the singular 
points will be very much abridged. 

The equation 5 = is an algebraic equation of the sixth degree ; the term 
with J as a factor is also of the sixth degree. The highest power of x, and also 
of y, is the fourth. The curve has the axes as double asymptotes, and has also 
two oblique asymptotes. The equation 5 = is of the seventh degree. As 
the highest power of x or oi y is the fourth, the curve has three asymptotes 
parallel to each of the axes. Among these are the two axes. There is one 
oblique asymptote. 

We come to consider certain restrictions which may be put upon the elements. 
In the planetary theory we are concerned usually with the case of orbits having 
small eccentricities. We proceed to apply this restriction to equations 5 = 
and 5 = 0, and to note what simplifications take place. It was found in the 
simpler cases that the abscissas of half of the singularities were of the order t or 
t', and that the other half were of the order 1/t or 1/t'. The same is true for 
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the ordinates of the singular points. If t, t' are sufficiently small we have 
approximately 

a- = 1 + T^ = 1 , a' = 1 + t" = 1 ; 

and if in addition x and y are small, we have approximately 

1-XT=1, (l-a;T)^=l, 1-2/t' = 1, {\-yT'y = \. 

In the expression for A the coefficients of (^y — t')^ and (1 —yT Y may be 
written respectively 

Y(X _ t)2 - y8(l - Xt)\ 7„(1 - XTf - /S„(X - Tf, 

which become, omitting terms of the second order, — /8 and 7^ respectively. 
Making in ^ ^ the approximations indicated and arranging in powers of y , 

we get 

P^^fl-^x{l+j)-\ + y-[{x-Tf + m,x'+'L^T'{\+j)x] 

+2/[-;S„a;(a;-T)^-2y8^„x^T'-;S/a,(l+j)+i7„a']+^y8„xV''=0. 
The equation i? = likewise takes a reduced form. We have approximately 

DA 8A ' 2 

^-^^ = %-^{y- ^')% ^ - 2/ :^ = To + ^(2/' - '^' )• 

Introducing these expressions in ^ = , simplifying as above, and arranging in 
powers of y , we get 

E, = y'[al3'{x - t) + c/3'cc] + y'[2cT(x - t) - 3c^Vx] 

+ 2/[— C(Sua;(a;^— t^) — a/8|,(a5 — t)^— 2cT'T(a; — t) — 2a^d^Tx{x— t)] 
— a/3'T' ( X — T ) + 3c/3't' x — Jc^^x +_/a7j(x — t) 
+ c/3ut'x(x^— t^) + 2a/3/3^ t' 'X ( X ~ t) — c/3't' x +_yc7gT'x= 0, 

where ,S' = ,S(1 +j). 

Thus we have a system of equations, (1), (2), P^ = , J?j = , for determining 
the singular points of <J>(») which lie near the origin. It is unnecessary to dis- 
cuss the geometry of the curves Pj = , J?j = . We know that in the region 
in which x, y are of the order t, t these represent very closely the curves 
P=0, P = 0. 

The equations for Pj , i?j are of the forms 

P^ = a(»2/« + aC^V^ + a(?J + af = 0, 
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where the coefficients d-P , ^p are polynomials in x in which j is the highest and 
i the lowest power involved. Eliminating y from these equations we obtain the 
determinant 



A, = 



4" 


af 
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This equation is of the fourteenth degree in a;, the lowest term being of the 
second degree. Thus a; = is a double root ; the other twelve roots give the 
singular points D, F, T, C, B^, B^, S^, S^, U^, U^, Ug, U^. 

The equation Aj = can be solved only when it is made numerical by giving 
particular values to the constants. The roots, which may be complex for general 
values of the coefficients, must first be separated and their values found by con- 
tinued approximations. This procedure gives, indeed, the singular points 
required but leads to no criterion of admissibility. A more fruitful method is 
to trace the original curves i^^ ^ , R^-=(i for certain values of the elements, 
and then, varying these elements, to note what changes take place in the singular 
points both as to their positions and as to their admissibility. 

The equation i? = is homogeneous in c, a; consequently a variation of 
the ratio cja may be applied very easily. But in practice the initial choice of 
cja will usually be sufficient. In fact, the important terms in the perturbative 
development are those corresponding to small values of 

m^v + m^^v', 

where v , v are the mean motions. We have approximately, for n very large. 

av -j- cv' ^= , 



from which 



a v ' ' 



The most important applications of the theory under discussion will be to the 
case of Jupiter and a small planet. The particular values of the elements are 
chosen in what follows with reference to this problem. Let L\^ sin <^, etc., be 
the elements as defined above of the orbit of one of the small planets and let 
L^, sin (/)', etc., be the elements of the orbit of Jupiter. Assume the value 
yS = 2, then the successive reductions of — a/c become 1, J, ^5^, etc. This choice 
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has the advantage of making all the singular points real in the case of zero incli- 
nation and coincident major axes, the necessary and sufficient condition being 
cja < . For the purpose of more accuracy in representing the curves ^ = , 
^ = by means of P^ = , ^^ = we shall assume t , t', j somewhat smaller 
than in the case in nature just suggested. We assume accordingly that 

a = -l, c = 2, |;8| = 2, t = .03, t'=.01, j=-.0005. 

If the perihelia coincide at the node, we have 

/3 = /3„ = 7 = 7„=2; 
if at 7r/2 from the node, 

/3 = /3o=-7=-7„ = 2. 

On introducing these values, the equations P^ = 0, i?j = become 

Pj = 2/'[-2a;] +2/'[5a;'-.02a;+ .0009] + 2/ [-2ar'+ .04a;2-. 0020a; rp-OOlO*] 

+ .0004a;- = . 

Pj = 2/'[ 2a; + .06 ] + 2/2 [ - .0036 ] + y [ _ 2x=' - .lOx^ + .0068a; ± .0030a; 

- .000024 =F .000030] + .04a;=' - .0008x'^ - .000016a; =f .000020a; = , 

the double sign arising through 7 •= 7^ = ± 2 . 

In the table on the following page are the equations resulting from particular 
values of x near the origin. The roots of the equation are multiplied by 1 00 and 
the corresponding coefficients and roots are set down. In the first column the 
values of x are given ; the two single and two double columns following give the 
coefficients of 2/^ y^, y and the constant term of the reduced equation. The three 
triple columns give the roots of the cubic in y. They represent respectively the 
roots for j = 0;j= -.0005, 7 = 7^ = 2 ; j = -.0005, 7 = 7^= -2. 

§4. The Criterion of Admissibility. 

It has been stated in § 1 that for a singular point of $ ( s ) to be admissible 
two singular points of F{z, t) must have come into coincidence from opposite 
sides of the contour \t\ = 1 . But F( a, t) is a multiform function of z, t oi 
an infinite number of determinations. For F(z, t) contains the expression P 
uniform mx,y , which are in turn multiform in » , < , as is seen from the relations 

z ^= x^ 6^ Ty"© *'' '' t =? x'^ a ''^ . 

The discussion for admissibility must then be made for each determination of 
F{z, t), that is, upon each sheet of the corresponding Riemann surface. 
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For the particular case t' = PoiNCARfi has overcome the difficulty. Using 
this value of t' in z above and solving for y, we get 

y = z" X " e '" ^ '^ ', 
in which it is seen that y is uniform in z^'" , x^'" , i. e., in s', x' if we put 

x" = x , z" = z . 

Also P{x, y) is uniform in z , x , and in these variables F[z, t) has only two 
determinations, namely those of the radical P~^ . This is effectively a reduction 
to one sheet, as these two determinations are equal in numerical value and oppo- 
site in sign. In the new variables the initial contour [ i | = 1 becomes jx' | = 1 . 
This method has been extended by Coculesco to the case of small values of t', 
where as a first approximation the eccentric and mean anomalies are considered 
equal. 

In general it can be shown by expanding z' as a power series in t', and neg- 
lecting terms above the Mth power (as is permissible except for values of y 
of the same order as t', such as occur along certain branches of the curve 
J^ ^ 0) that F{^ z, t) becomes a function of multiplicity 2n in the variables 
z', x . The geometrical meaning of this transformation as applied to the con- 
tour 1^1=1 is interesting, but the investigation is unnecessary to the present 
discussion. For it will be shown that the admissibility of points for the general 
case of t' =1= can be determined from the discussion of the case of t' = . 

In fact, PoiNCARfi has shown that, in general, a singular point of ^{z) can- 
not change its character as regards admissibility for any variation of the ele- 
ments. The exceptional cases, as we shall see, are those for which there arises a 
change in the combinations of the singularities of i^(K, i5) which form the singu- 
larities of <I> ( a ) . We proceed to establish a theorem simpler in form than the 
original theorem by Poincab:^, but embodying its essential point. This theorem 
may be enunciated thus : 

I. For the values |a| = 1 the singular points of F(^z, t) cannot cross the 
circle \t\ = 1 . 

Having seen how the circle | ^ | = 1 can be replaced by either | cc | = 1 or 
|x' I = 1, we are at liberty to make the discussion in the plane of x. Consider 
the relation 

z = x^e''^-' ^ ^ y'e"' ^'J "'. 

If we let 2 be a parameter, this relation is an implicit multiform function of the 
two complex variables x and y. We have seen that it has an infinite number 
of determinations. There corresponds then to the function z(x, y) a Eiemann 
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surface S with an infinite number of sheets. The surface varies when the par- 
ameter z is varied, and for | » | = 1 becomes the surface S^ . Upon this surface 
trace the circle C^ defined by 

l^'l = \y\ = 1- 

This is possible, for if x be given a value of modulus unity then we can choose 
a value of y having the modulus unity in such a manner that z may have any 
value whatever of modulus unity. On the other hand if | x | = | y | = 1 and 
X = y, then there corresponds a unique point on the circle of unit radius in the 
g-plane. 

Now to each point on the Eiemann surface S there correspond definite values 
of X and y , and consequently definite values of the eccentric anomalies u and u. 
The points on the circles 

give real values of the eccentric anomalies, i. e., a real point on each of the 
orbits. But in this case F(^z, t) has a singularity only when the distance 
between the two bodies is zero, that is, when 

or, when the distance of one of the planets from the sun is zero, that is, when 

1 — sin (^ cos M = , 
or 

1 — sin <^' cos u' = . 

Excluding then the case of intersecting orbits and also the limiting case of the 
problem of two bodies where the orbit reduces to a straight line, it is impossible 
that for any variation of the elements a point on the circle C^ shall become a 
singular point of F{z, t). Therefore a singular point of F[z, t) can not cross 
the circle | i | = 1 for values \z\ = 1. 

A singular point of ^{z) has been defined as being admissible when it is 
formed by the coming into coincidence of two singularities of F{z, t) from the 
interior and exterior respectively of the initial contour | i | = 1 . The theorem 
just proved may be applied to the question of the permanence in the character 
of admissibility of the singular points of 3> ( a ) . To distinguish among these 
singular points certain definitions are useful. We assume that the elements of 
the orbits, with the exception of the anomalies, remain constant. Let A , B 
be two singular points of ^(s) at which z' has the values z[, z'.^ and assume 
\z\\ < I ^2 1 < 1. When z is made to vary starting from the point A , then the 
corresponding points in the plane of x separate and move subject to the condi- 
tions defining singular points of jP(», t). If in the plane of x there is a 
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means of direct passage from A to S along a line of continuously increasing 
modulus of z, then B is subordinate to A.* The means of passage is the locus 
furnished by one of the equations (1), (2), P = . 

If the point z encounters no other singular points in its course, and if the 
points A and B are only double points, then A gives three terminal values in 
the plane of x, one proceeding from A directly and two proceeding from the 
branch separating at B. 

We define a singular point as primary if it has no admissible subordinate 
point. It is seco7idary if it has an admissible subordinate point. A primary 
point may be either admissible or inadmissible ; a secondary point is pseudo- 
admissible, for, although there correspond to it terminal points both within and 
without the circle |a;'| = 1, there can be no change of elements which will bring 
it upon the circle of convergence \z'\ = r without also causing it to pass through 
the circle of absolute value of its subordinate point, and the subordinate point 
itself determines the limit of convergence of the Laurent series. 

With these definitions we proceed to the discussion of admissibility, taking 
up the different cases that may arise. It is requisite first to establish the 
theorem : 

II. A primary singular point q/" <!>(») can not change its character as 
regards admissibility for any variation of the elements, and this character can 
be determined by varying z along any path whatever of increasing modulus. 

Consider the method of formation of the singular points of ^{z). The vari- 
ables X and y satisfy certain equations of the series (1), (2), i-* = 0, i? = 0. 
The intersections furnished by these equations give singular points of <i>(s) of 
the first and second species, and to such a point there corresponds a definite 
value of z\ that is, a point in the plane of z. Instead of finding the intersec- 
tions indicated and then computing the values of a', it is simpler to introduce 
certain auxiliary curves to take the place of (1), (2), P = 0, M = Q, whose 
intersections give the required singular points of •^'(a). Let 

a{x, y) = 

be any one of the curves (1), (2), P = , i? = . We have 

When r varies in a given manner, the points of intersection of a = and z = / 
describe certain curves ; in particular, if z describes the circle 



the singular points arising from the intersections of the curves will cause x to 
*POINCAEE, loc. cit., p. 308. 
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describe certain /3-curves in its plane. The curves a = may be replaced by 
the corresponding /3-curves. Mere intersections of the /3-curves are not singular 
points of 4>(a); the tracing points must coincide for some value of a'. 

For the initial circle | a' | = 1 the /3-curves in the plane of x do not intersect 
the circle |q3'| = 1, by theorem 1. Now let the circle |a'| = r' contract, keeping 
the origin fixed. The /3-curves sweep over certain portions of their plane. The 
significance of this is that corresponding to any continuous path whatever from 
any point on the outer circle | a' | ^ 1 to any point on the inner circle | a' | = r' , 
there are certain paths lying entirely within the region swept over by the /3- 
curves. Suppose that when the circle | a' | ^ r' is reached two tracing points of 
the /3-curves coincide. Corresponding to this point in the x'-plane is a point on 
the circle | a' | = / , namely one of the singular points of <E> ( a ) . Now reverse 
the variation of a', and starting with the singular point let a' go to a value 
I a' I = 1 along any path whatever consistent with the continuous increase of the 
radius of the circle on which it is found. The points in the oj'-plane accordingly 
trace definite loci in the swept-over area and rest finally on the initial /3-curves 
corresponding to | a' | = 1 . These curves being entirely within or entirely with- 
out the circle \x\ = 1 , the terminal values of x are definitely within or without 
the circle | x' | = 1 for any a'-path to the circle | a' | = 1 ; whence we conclude 
that the point can not change its character as regards admissibility for any 
change of the elements, for this would necessitate that some of the points of the 
/3-curve should cross the contour 

which we have shown to be impossible. 

If in a series of singular points, connected by relations of subordination as 
defined previously, the one first encountered is inadmissible, the circle | a' | = r' 
may be contracted over this singular point, since a deformation of the contour 
in the plane of x is not rendered impossible. This may be repeated until an 
admissible point is reached. It is evident that the theorem applies to this point 
as well as to all the preceding points of the series, that is, to all the primary 
points of the problem. 

Consider now all the points both primary and secondary in the a-plane con- 
nected in any manner whatever as required by the definition of subordination 
and as permitted by the method of varying a . If the elements be varied con- 
tinuously these singular points likewise vary continuously. Their curves of con- 
nection must be deformed if necessary to correspond to the displacements. It 
is proposed to account for the behavior and the consequent character of the 
points for all possible changes on their relative positions. 

Let A represent a typical singular point with B as an immediate subordinate 
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point, and consider any particular series of connected points. Let B pass 
through A . Then : 

(1) If both were originally inadmissible they remain so after the change, for 
their combined terminations in the plane are still on the same side of the unit 
circle. 

(2) If A were originally admissible and B inadmissible, the passage of B 
through A leaves A admissible and renders B admissible and secondary, i. e., 
pseudo-admissible. 

(3) There is no case of B admissible and A inadmissible. 

(4) If B is admissible, then A is pseudo-admissible, and the passage of B 
through A gives the following result : 

(a) If A were intrinsically admissible, that is, if its admissibility did not 
depend upon a branch of B not common to A and B, then A remains admis- 
sible and B becomes pseudo-admissible. 

(6) If A were not intrinsically admissible, then A ceases to be admissible, 
and B remains admissible. 

This discussion covers all cases of interaction of primary singular points. 
This interaction produces certain secondary points. If we consider the passage 
of a secondary admissibility through an immediate subordinate primary admis- 
sibility, we shall have one of the above four cases, consequently the converse 
results hold. In fact, these converse relations are contained in the cases as 
already stated. We have : (1) is its own converse, (4) (a) is its own converse, 
(3) does not exist, (2) is the converse of (4) (6), 

It has been assumed that B actually passed through A, but the results would 
be the same if it passed by it simply, i. e., passed through it in absolute value. 

It should be remarked here that the variation of z from a secondary point to 
to I s I = I is not entirely arbitrary. For if all subordinate points are not 
taken account of, it may easily happen that we get an apparent inadmissibility 
when the point is actually admissible as found by the path through the sub- 
ordinate point. 

This completes the theory of admissibility. Before proceeding to the discussion 
of the points found in § 3 it should be remarked that there are for each point 
in addition c — 1 others of the same modulus and differing in amplitude by 
the arcs kirjc. Poincar:^ showed that these points are all admissible or inad- 
missible according as one of their number is admissible or inadmissible, and 
this holds true in the general case. 

§ 5. Discussion of the Singular Points. 
Having seen that the singular points which arise in case of particular elements 
can not, as a rule, change their character as regards admissibility for a change 
of elements, it is advantageous to make use of the simpler cases to determine 
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what points are admissible and to follow these points when certain of the ele- 
ments are made to vary continuously, in order to determine in the general case 
the points that are admissible. 

Accordingly we make direct use of the results of Poincar^ and Coculesco 
already cited. Let us designate three cases as follows : Case 1, j = ; Case 
II, J =1= , perihelia coincide at a node ; Case III, ^' 4= , perihelia coincide at 
ninety degrees from a node. In the figures 1, 2, 3 the light and heavy lines 
are the curves in the case J = , J 4= ^ respectively, also the full lines are the 
curves 1*^= 0, x= t, y = t, and the broken lines are curves iJ^ = 0. The 




Fig. 1. 



positions of the small arcs show where, with respect to the circle | a;' | = 1 , the 
points on the curves Pj = , x = t , 2/ = t' lie for which \z'\ = 1 . 

Case I. J = 0, « = 5'. This is the case discussed by Poincar^ and 
Coculesco. The singular points have already been classified in § 2. Of the 
points of the first species B , S, £J, W, X, Y are inadmissible and U, V are 
admissible. Of the points of the second species, C, F are inadmissible and Z>, 
:7'are admissible. 

The variation of z is as follows: On the branches ^^a^■, z decreases to C, 
which in this case coincides with 5 , is a minimum at C, increases from C to O , 
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changes at O from |2!| = ooto|a| = 0, and increases from O through E. On 
the branches /S^/Sg, 0L^a^,z decreases to Z), is a minimum at D, increases through 
B toP, |s| = oo at P , |s| = Oat Q, increases to F, which in this case coin- 
cides with i?, is a maximum at F, decreases from F through V. On the 
branch ^^a.^, z increases to T, is a maximum at T, decreases from T to O, 
changes at O from |»| =0 to |is|^ oo, decreases from O to U, and increases 
from U. Along a; = t the value of z increases with y positive, and along y^r 
it increases with x positive. 

Case II. j' ^= , w = ftj' = . The curves in this case are given in Figures 
1, 2, which are plotted with the values given in Table I. The branches a^, a^. 
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Fig. 2. 

etc., are the curves of the Case I, just considered ; the branches a[, a'^, etc., are 
for present consideration. The curves Pj=0,x = t,2/=t' are represented 
in full lines, and the curves i?j = by dotted lines. 

Consider the changes which the original curves undergo when the inclination 
is made to vary from zero, in other words when j varies from zero to a small 
negative quantity. The segments a^ R and It O vary continuously to form the 
branch a[0. The segments a^R , RU, Ua^ form the branch a'^U^ctl. The 
segments O U, Ua^ form the branch OU^a.'^. The straight lines a? = t , 3/ = t' 
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are of ■ course unchanged. The branches a^ a^ and a^ a^ become the branches 
ajOtg and oL^a'^ respectively. 

The singular points for y sufficiently small will lie very near those of Case I. 
The point R is in reality a combination of what would be in general the points 
F and li . This coincidence is due to the choice of constants in the particular 
numerical case computed. A change of the ratio cja is sufficient to separate 
the points F and i? without changing the curve ^ = on which they lie. But 
these points can be separated otherwise. The branch a^RU is tangent to the 
branch a^RU; this gives two coincident points of intersection at i?. The 
branch a, i?Z7a]so cuts the branch a^iJO at i? in a point of the second species 
R^. Thus R is composed of the points F, R^, R^, all of the second species. 
When j varies in the manner indicated the point F remains real and takes a 
position F^ ; the two points R^ , R^ become conjugate imaginaries. 

The point U in the simple case is composed of nine points. These arise as 
follows : each of the two branches oi R = which pass through fj cuts each of 
the branches of P = at the same point, giving rise to four points fjj, U^, U^, Cf 
of the second species ; each of the straight lines cc = t , y = t cuts each of the 
branches of i-* = in the same point, giving rise to four coincident points of the 
first species, t^ , C^ , Uj, U^; the point U is also the intersection of x = t, 
y = t' . Now vary j as before ; the nine points all separate. U^ and C^ of the 
second species remain real. U^ remains real. U^ is imaginary in the figure, 
but it is evident that for sufficiently small values of j it is real. In fact it 
leaves C^ as a real point and different from the point E, these two points 
approaching each other along the line a; = t . For a certain value of j the 
points F and U^ coincide, and when j is given still larger values they become 
conjugate imaginaries. The points U^, U^ are conjugate imaginaries; the 
points C^, C^ are also conjugate imaginaries. 

The remaining points F, of the figure, and W, T', X, not represented in the 
figure, take positions near their respective original positions. 

Discussion. — By virtue of the theorem established in § 4 it is not necessary 
to discuss the points E, F, R^, R^, V, W, X, Y, since we have already 
determined the nature of the corresponding points in the simpler case. We 
have only to consider the nine points corresponding to U. 

The equation /* = is of the form 

P=ST + jA. 

We proceed to find the form of this curve for the values of x lying between 
some value x of the order t , and the value a; = 1 . The product 8T considered 
as a polynomial in y is of the fourth degree, while ^ is a polynomial of the 
third degree in y. By taking _;' sufficiently small we can make the curve 
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lie as close as we choose to the curve 

ST=0 

throughout the region under consideration ; for the coefficient of the highest 
power in y remains unchanged, and all the coefficients of A are either finite or 
small. 

Then for a given value of x the corresponding values of z along similar 
branches in the two cases will be approximately equal. In the case of j* = 
the terminal values of x, when z described the straight line of constant ampli- 
tude, were definitely within or without the circle | a;' | = 1 . It is possible in the 
present case to choose j small enough to make the terminal value of x', relative 
to the circle | as' | = 1 , correspond exactly to those in Case I. The character of 
the real points Cj , W^, U^, U^ is now evident. 

f7g is formed by the intersection of a; = t with a'^ U^ a'^ . The two points on 
separation follow their respective branches in the direction of y increasing. They 
reunite in £1. But £! is inadmissible, therefore U^ is inadmissible. 

C^ is inadmissible. The two points which separate are found to follow 
fT'j O and a; = t in the direction of y decreasing. The terminal values of x are 
both without the circle | cc' | = 1 . 

Z7j gives z' a minimum. When z' varies along the positive part of the real 
axis of its plane toward the value z' = 1, the two points x' on separation move 
along U^O and C^jK^ respectively. The latter encounters in its course the 
point U^ . All the terminations for Z7, have been found to lie within | a;' | = 1 . 
The same is true for the terminations on the branch Z7, a^, therefore ?7, is 
inadmissible. 

ZZj is admissible. Since at this point z' is a minimum the values of x' on 
separation follow the branches U^a^ and U^o^'^. When z' attains the value 
unity we have for the first point x' < 1 and for the second a;' > 1 . 

The points U^, U^, U^, U^ are left for consideration. These points are aU 
imaginary but their amplitudes are very small when j is put small. Correspond- 
ing to each of these points the value of z is imaginary with a small amplitude. 
Instead of varying z' along a line of constant amplitude, as was done in the case 
of the real points, let us vary z' along the straight line limited by the points z'^ 
and »' = 1 , where »^ is the value of z' at one of the singular points. The loci 
of the singular points which separate when z' is made to vary from the value z'^ 
are certain curves satisfying the relations 

for points of the second species W^, U^; or 

z{x,y) = ^, y — t'=0 

Trans. Am. Math. Sec. S3 



504 F. w. reed: on singular points in the approximate [October 

for points of the second species, U^, U^, where ?" is a parameter varying along 
the straight line limited by z'^ and «' = 1 . In the case of J = this variation 
of f is along the real axis, and consequently the singular points on separation 
follow P{x, y) = Q m the real plane. For j very small the singular points will 
describe curves lying not in the real plane, but very near the real part of J' = . 
The approximation becomes closer as ^ approaches f = 1 , for which value the 
singular points have positions on the real part of P = not far from those 
terminal positions found for j = 0. 

Consider the singular points of the second species, U^^U^, formed by the 
intersections of the curves P = and J? = in the neighborhood of the 
point U. Such an intersection implies that z is a maximum or a minimum. 
In each case we know it to be a minimum. Then the four points x after 
separating, two from U^ and two from U^ , follow each of them one of the four 
segments of P = , which approach coincidence with UO , US , Ua^ , Ua^ , 
respectively. 

An examination of the order of interchange of the branches near U is neces- 
sary. But this discussion is practically obviated through a priori considera- 
tions. In fact, U represents an ordinary node of the function P = in the 
case of j = 0. When j varies the interchange may take place in one of four 
different ways. The double point may continue to exist near U, or the two 
branches may cease to intersect and the continuous branches may lie near the 
original curve in any of three ways as follows : a'^ UR and a'^' UO, a'^ Ua'^' and 
R UO, a,' UO and a" UR . It can easily be shown that if the node continued 
to exist for variations of j, there would be one or three admissible singular 
points according as one of the two points of the second species is or is not 
on the segment U^^^a'^. But it can be shown that for general values of j 
there is no node C/j^. Consider the three other possibilities. By virtue of the 
continuity of the function P we see that each of the four segments which 
the singular points describe after separation is involved once and only once. 
We know that of the four terminal values of x three are within the circle 
I a;' I = 1 and one (on the branch a'^) is without this circle. Thus of the two 
conjugate imaginary points, U^, U^, one is admissible and one is inadmissible. 
This distinction will be considered later. 

The points U^, U^ which are the imaginary intersections of y = t and P = 
in the neighborhood of U become real when <!) = «' = 7r/2 , and are reserved 
for discussion under Case III. 

The point U, which is the intersection of a; = t and y = t', is or is not 
admissible according as one of the points U^, Ug, supposed subordinate, is or is 
not admissible. This illustrates a distinction made in the foregoing section. 
In fact, if the points on separation were permitted to follow the line y = t only 
in the real plane, we should be led to conclude that U is necessarily inadmissible. 
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The curves as they exist in the third quadrant are represented by Fig. 2. 
The branches /3j, /S^, etc., are the curves of Case I, while ^[, 13'^, etc., are the 
curves in the present case. From the figure it is seen that as j varies from zero, 
the line fi^TOB^^ becomes ^[^'^, and /3^DB^^ becomes yS^/S^. The branch 
0^DTOB^^ becomes ySgO/S^. At B three points coincide. These separate 
when _;' is made to vary, the point C remaining real, B^ and B^ becoming con- 
jugate imaginaries. The points D, T, remain real and their abscissas increase 
in absolute value. 




The points B^, B^, C are inadmissible and the points £>, T are admissible, 
as was found in Case I. 

Case III. J 4= 0, « = 5' = IT 12. In this case the curves as they exist in 
the first quadrant are represented by Fig. 3. They are plotted with the partic- 
ular values given in Table I. The branches a^, a^, etc., are the curves in the 
case J = 0; the branches a[, a'^, etc., are for present consideration. 

When j varies from zero to the small value used in the table, the original 
curves vary as follows: a^RcL^ becomes a'^a'^; ORUO remains a closed curve 
and recedes toward the point O; o.JJa^ becomes a.'^o^\ the branches a^RUa^ 
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and a^Ua^ become respectively the branches a[RUa'^ and a'^a'^; finally, the 
branch a^ passes through a loop form, then to the position Oa'^ . The order 
of the separation of the branches here is seen to be opposite to that of Case II. 

The singular point at li breaks up into the points F, S^, R^. The points 
F, R^ are conjugate imaginaries. R^, the point lying below R, is real at first, 
but for the value of j employed in the figure it has become imaginary. When 
the loop formed by the segments OR UO begins to recede toward the origin, 
the two points of intersection of this loop with a line R U, little different from 
the lower of the dotted lines RU, will be real and different. Designate by U^ 
that one of the points near U. The value of j continuing its variation, the 
points C/j and R^ come into coincidence and then separate to become conjugate 
imaginaries. Consider the variation of »' on this loop when U^ and R^ are 
real. The loop is cut by the line y = t' at the points V and U^ . The varia- 
tion of »' is as follows : zero at O ; O VR^ , increases ; maximum, at R^ ; 
R^U^, decreases; minimum at Z7^ ; U.^U^O, increases; infinity at O. Now 
let the loop decrease, the maximum point R^ comes to coincide with the mini- 
mum point C/j, then z' continuously increases around the loop OV^U^O. The 
loop continuing to decrease, V^ and U^ approach each other, coincide, and then 
become conjugate imaginaries. The curve a[ TJo.'^ and the line y ■= t' cut the 
curve J° = O at approximately the same point in the figure. The straight line 
a? = T cuts the curve -P = O in two imaginary points U^, U^. 

Discussion. — By using the results of the two preceding cases the character 
of the singular points in the present case is easily determined. 

F is inadmissible. 

F is inadmissible. 

R^ has F and F as subordinate points and is itself inadmissible. 

R^ is a maximum point of »', and when z increases the two points become 
conjugate imaginaries at first and then coincide at R^. Thus R^ is inadmissible. 

U^, the intersection of a'^a'^ and Ua'^, is admissible. For, after separation, 
the two values of x follow the branches a'^ and a'^ respectively, the terminal on 
the former branch being without and that on the latter branch being within the 
circle \x' \ = 1 . 

f7^, subordinate to U^ for sufficiently small values of J, is admissible. The 
tangent to the curve Ua'^ at U is negative thus giving V^ a position above U^. 
The terminal value oi x in the directions of x increasing and a'^ are respectively 
within and without the circles | a;' | = 1 . 

U has C^ as a subordinate point and is consequently admissible. 

U^ is admissible ; of its two branches one follows the curve U^ , O and the 
other passes through U and C^ which are admissible. 

U^ is admissible ; it has R^ and U^ as subordinate points, of which the latter 
is admissible. 
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Fj is admissible ; it has at first M^ and U^ as subordinate points, and later U^ 
alone. 

U^ and U^ , the intersections oi x=-r with P = O are admissible and inad- 
missible respectively. For, of the four branches UR, UO, Ua^, U<x^, only the 
latter gives a terminal value of x' greater than unity. With the argument used 
for CI and U, of Case II, the conclusion is evident. The discussion for the 

3 4 ' 

third quadrant is not difficult and may be omitted. 

§ 6. The Case of t' = , J 4= . Conclusion. 

Before passing to conclusions relative to the general case, it is interesting to 
consider what simplifications arise on taking the eccentricity of the orbit of the 
second planet zero, the eccentricity of the orbit of the first being small, and the 
mutual inclination small. 

Neglecting in the equation JP = the terms in t' and arranging in powers of 
y , we have 

y^lyajx{x-ry-^,7x{l-xrf{l+j)]+yl(x-Ty(l-XTf+^^yx'] 

+ y,ojx{l - xrf - ^,ax{x - t)^(1 -fj)] = 0. 

An approximate expression for this curve in the neighborhood of the origin can 
be obtained in a very much reduced form by neglecting all but terms of rela^ 
tively low order in the various coefficients. Thus considering x of the same 
order as t, and neglecting all but terms of the second order in the coefficients, 
there results the equation 

P, = f[-0ax{l +i)] + y[ix-rf+ ^^„<rW] + %ojx = . 

Likewise the equation H = becomes 

i?, = f[al3{x - t)(1 +i) + c0ax{l +j)] + y [2cT(a!- t)] 

+ [-a^o(»'-'^)'-<'^o'^(«''-'^')+i{- C7oO-»'+ a7o(a'-T)}] =0. 

The eliminant Aj = of Pj = and H^ = may be found as in § 3. It has for 
its roots the abscissas of points corresponding to U, F, R^, R^,D, 0,B^, B^. 
Equations P^ = , R^ = may be written in the expanded forms 

^2^ y'[«ll»'] + y [«22"'' + ^21^ + «2o] + [»33«'' + h^^^ + hl^'\ = ^' 
i?2 = y'^ihl^ + «lo] + y ihl^ + «2o] + [^33^^ + ^32"'' + t'zl^ + C] = 0, 

where Sjj, t'^^, t'^^ are so designated because they contain all the terms in j, and 
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to these terms the double sign is attached. Then A^ becomes 

[Sjja;' + 83233^ + SjjCb] [MjCc' + u*x* + • ■ • 4- m,] 

+ [tsB^ + ht^ + «> + «3o] 1%^' +■■■+'»,]= 0, 

where the u , v are functions of the s , ^ as indicated. Or again, 

Aj = XgX* + \x' H H X„ = . 

Assuming the approximate values of the roots as found by plotting, they may be 
corrected by applying the well known formula 

which is true for all simple roots, real or imaginary. The work may be reduced 
in the two cases under consideration by finding the four real roots first and then 
reducing the order of the equations. 

In the two cases of w = m = and 5 = &>' = ■irl2 , the roots of A^ may be 
calculated very rapidly by the above formula when j is varied little by little, 
since changes in _; will affect Sjj , ^3, , t'^f, and u^ , v^ in an additive manner simply. 

Since the choice of /3, t can be made once for all, and since j can be made to 
assume any assigned value in the manner indicated, there remains for considera- 
tion the variation of the longitudes of the perihelia. Having taken t' = , we 
have a> = a and 

L'^ L'' - 

/3 = /3„ = ^-,, 7„ = ;pe-*"- 

When a is taken different from zero, imaginary terms affect the expressions X. . 
The roots of A^ = become imaginary. They may be computed by the above 
formula as m varies from zero. The computation of af^ in the functions 
Aj(Xd), ^\x^) is best effected by the use of trigonometric tables by writing 

fljjj = |a!|,|"e*"* = |a;(,|''(cos n6 + i sin n9). 

From the foregoing case the means of extension to the most general case is 
now evident. We have seen that it is necessary to take account of the analytic 
possibility that a singular point may become admissible or cease to be admissible, 
but slight attention to the geometry of the variation of the singular points as 
the elements vary will suffice to guard against error. Probably the selection of 
the primary admissible points, as well as their tracing in the x-plane, with the 
application of the principle of the limits of the roots of an equation, will answer 
in most of the practical applications. Approximations may usually be made in 
the manner employed here, and where necessary recourse may be had to the 
original equations of a higher degree. 
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Only points within the circle | s | = 1 have been considered in the foregoing 
discussion. After we have determined the coordinates x,y oi the singular point 
which gives z nearest the circle, the point outside the circle | as | = 1 and nearest 
it is found as follows. We have 



A(jc,t) = XA(^^, -j) = 0, 



where ^ is a power of x and A is the equation in x determining singular points. 
This shows that \jx is also a root of A if — i is put for i. But this is effec- 
tively a change of elements, namely an interchange of /3 with ^^ and of 7 with y^,. 
Make this change, then 



A(x,-i) = XAQ, i) = 0. 



Thus we see that by putting — i for i the knowledge of the inner point which 
conforms to the demands of the problem enables us to find the required outer 
point. 

The simplest relations from which to begin the variation of the perihelia are 

(0 = 0, 2 5 1", -^, 

with 5 coincident or opposite, which gives 

^(=/3„) = ±7(=7o) = ±X^- 

The variation of c and a has not been considered. These constants enter the 
equation H = only. This equation is homogeneous in c , a, so that if it is 
decided to vary the ratio c/a after a singular point is found, a factor may be 
entered readily in each term. Otherwise the problem may be begun anew for 
each choice of c/a, c and a being small integers. 

The Univeesity of Illinois, 
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